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ABSTRACT 


A continuous identification of parameters is performed 
on a simulated fast breeder nuclear reactor system using 
hybrid computation and applying techniques of statistical 
regression analysis and exponentially-mapped-past functions. 
Output states which are not directly measurable are estimated 
by use of a Kalman filter. 

The method developed in this study is applied to a nu- 
merical example which demonstrates that unknown parameters 
can be identified within 3% of their actual value, with sig- 

r 

nal noise ratios as low as 10:1 in the measured states. The 
example also demonstrates that convergence occurs in a rea- 
sonably short time. 

~ An executive software routine was written in FORTRAN IV 
mer the XDS-9300 digital computer and was applied both for 
mies control of the process and tor establishing the opera- 
MPiomal performance and accuracy of the Comcor Ci-5000 analog 


Computer that simulated the nuclear reactor system. 
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I. INTRODUCTION 


at 


Fast breeder nuclear reaetons such tas the srecenel, 
proposed design described in [1] do not have inherent load 
following ability and therefore require an automatic con- 
Mmeorier to handle load changes. It is important in désign- 
ing for stable operation of the automatic control system to 
know the value of the parameters of the plant within a rea- 
sonable degree of accuracy. These parameters change due to 
such factors as temperature, fuel depletion, fuel loading 
ead power level, [2], making it desirable to have continuous 
identification of parameters. 

The continuous identification of parameters of dynamic 
systems has received intensive attention during the last 
mew years [3], [4]. References [5] and [6] among others, 
have shown successful results by applying gradient techniques 
ama quasilinearization, respectively. These methods are cum- 
bersome for nonlinear systems and in particular for dealing 
feeene a nonlinear space-dependent reactor system, [6], where 
Meemrequirement of continuous identification cannot be met 
eamice the in-core sensors which are necessary for certain 
state measurements have not yet been developed’. 

The implementation for continuous parameter identifica- 


tion is difficult and complex. For example, many gradient 





1 

The procedure available consists of the exposure of 
a matallic foil by immersion and takes enough time (between 
one and two hours) to eliminate any possibility for con- 
UEMUOUS Identification Of parameters - 








techniques, even when they represent in some cases very ele- 
gant ways of approaching the hill-climbing search procedure 
are limited by practical considerations in theimeanplaeaetene 
lo 2zilustrate this, reference 15 made fo a partreula mone 
ent (steepest descent) technique which uses an error crite— 
rion and sensitivity equations, and which, among other 
mamgerrements, has the need for a forcing function. A block 
diagram of the algorithm is shown in Figure 1. In this 
method it is required to derive, implement and solve the 
Semsitivity equations in order to obtain the value of the 
meaarent for the specified cost function. In order to 1den- 
tify n parameters the number of analog components, using 
Seandard methods, must be multiplied by at least n * | Samnce 
mowed linear case these equations take the same form oi the 
moet, [7]. In addition, for nonlinear systems the Sensi- 
tivity equations have more complicated forms which do not 
mec the use of techniques such as mul€iplexing, [7]. To 
apply these methods the system must be initialized with some 
maev tous information, usually an educated guess of the pa- 
fame ter Values, and then the process must be iteratively 
meme tialized until some criterion 1S Satisfied. Since the 
parameter correction is a function of the gradient and 

Smee Chis gradient approaches zero as the parameter value 
approaches the actual value (solution), the procedure will 


converge in a reasonable time only if a means is provided 


More recent techniques, [8], allow for simplarresuie 
when the model can be expressed in canonical form. 


EO 
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to speed up the convergence sin ties tei) sot seme iianiay 
Ope the cost function. 

Other techniques, [9] and [10], which are used in non- 
Mimear systems suffer fromemore or less the Sanem@emeb lence 
Mimese difficulties ondieate the need fomean 1demer preaunen 
eyoetem that will operation without supervision, require 
mearencrmcuuecatcd@ptesses Sor eini tmalwzation, be very reli- 
able, and use simple implementation that avoids the trouble- 
Meme iterative initializing process. 

The procedure developed in this thesis makes use of 
mme method of statistical regression analysis, applied pre- 
viously with success to linear systems, [4], and confined 
ter a long time to the analysis of statistical populations, 
mee) and [12|. The method of steepest descent is applied 
memeenc classical regression analysis in order to speed up 
the convergence. To avoid the need for reinitialization, 
meee properties of exponentially-mapped=past statistical 
Pemeaables 1s used, [13}. 

MiENOigheene DhOGcedunre 15 japnlied Nere Slama ss Imp lsmimled 
mocel of a fast breeder nuclear reactor system, the method 
memeeneral enough to be applicable to many other automatic 
Bemtrol systems that are represented by more complex models. 

Forethie particular muclear reactor system used as the 
example application, the variable representing the delayed 
mewcron Concentration 1S a nonmeasurable state and in order 
Hemmave ItS value for the identification system a Kalman 


Peter sl) ana /15], was used. 


I 








The unknown» parametersywerecetdemerii1cd Wore neo cece 
their actual values. The example demonstrates that con- 
Mereence occurs and that fhe time ot convergence conpane— 
mayorably with current methods, [16]. @@ihe simmikaeienmne- 
sults shown in this thesis were obtained using zero as the 
initial values for the parameters; however, tests have 
Bwewn that the convergence time can be greatly improved iff 
eamcducated guess is accepted. 

ATeexeGulive SOLtWaLe DTOWEINGE was WEleechurn eon n, 
meecor the XDS-9300 digital computer and was applied both 
Momethe control of the process and for establishing the 
Operational performance and accuracy of the Comcor-5000 
analog computer that simulated the nuclear reactor system 
ela was used to implement the analog part of the identifi- 


@ea1on procedure. 


ES 








II. REACTOR DYNAMICS 


The point model of thewniticlean neaeror kine cicecaud. 


235 
mons for U fuel can be found in several references, 


BiG] - [19]. 
These equations are 
6 

ee Ke 3 

Me ae ee ee ae (1) 
and 

os 

at = ren - AGS, de Ae Oe (CZ) 


The left-hand term in equation (1) represents the rate 
Seemeiiange of the neutron density; the first term on the 
mee Side of equation (1) represents the rate of production 
Paeprompt neutrons, the second term the rate of absortion 
plus leakage and the last one, the production rate of de- 
layed neutrons. The variable 6K, the reactivity input, can 


bemexpressed as 
i 
OK = Oe >. Oeless. 
qe Jog 


Applying this substitution, equations (1) and (2) become 


QR 6 
dn _ - 3 els » 
ar =|. ra asTj | qe gent & Ace (3) 
oT ee 
aie ge ce Oy, 
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where T is defined by 


. dT. (5) 
| n 


be ut 2 - a] - yee it, at i=1,...,2. 
= t Sy Ny | k=] jk | J «| ae 


Samce Us and 8, are fractional quantities. 1t follows tmem 


their definitions that 


Y 
yu, 71 
j=. 3 


ee 


1=] 


and 


mieminitial conditions for equations (3) thru (5) are 


It 
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reo) 





It 
®) 
i 
ee 

O 
i 
i 
-— 
Or 


C; (0) 


It 
am) 
iH) 
pee 

Ww 
Pa) 


T. (0 
; 69) J 
Hom the purpeee Gmthis Study a single-time-conseant 


meedback reactor with a single group of delayed neutron 


precursors will be used; the following equations define the 


Sys Cell 
gt > [oe - ot) gP- pha + 2c 
Se = sen - Ac (7) 
ec - =a E: - 1 anlar (8) 


1S 








ine initial’ condgditions -2nc 


(0) = 0 
n(0) = ny 

B Ny 
CO) = Co Tah 
T(0) = 0. 


By changing variables as indicated in [16], (see Ap- 


pendix A) a more suitable representation is 


N = 7& [(K, - aT - 1) N+¢ DI (9) 
e 
D=. (N - D) (10) 
T= yee Tomas (11) 
where 
Me = pent x = me Baca OL qg = & =e 
ie 2 Oe B 7 gk 9 re 
) 
eae 


mmo the initial conditions are 


N(0) 


iI 
-_ 


T(0) 


K (0) 


NH 
a) 


D(0) 


il 
bt 
ll 
= 


Experience with present day design of nuclear reactors 
indicates that these systems are stable under normal opera- 
mine COnditions. The large breeder reactors are still in 
the design stage and a computer simulation analysis has 
Bnown that these plants will be stable, [1]. Nonetheless, 
a stability analysis applying the method of Liapunov has 


been performed for this model and is included in Appendix B. 


16 








The model defined by equations (9) thru (11) was im- 
plemented, but acceptable accuracy could not be obtained. 
The large value of By¥e* required’ a Mamce, scaling -sactormeanG 
also a large gain for the state variable N. Because of 
mmese two reasons the voltage that represented the state 
variable N become too small and too close to the uncertainty 
Mmrel (fourth digit) of the analog computer. 

ErOm, Simulation te tse opsenveda tmat N 1S approximately 
pamlal to zero. This and the fact that 8/2* is very large 
gmplies that equation (12) is a satisfactory approximation 
meme equation (9). This further simplification will, of 
course, introduce some inaccuracies since it 1S assumed, in 
fact, that N is zero; however, it has been shown, [16], that 
Mies largest error introduced with this assumption is less 
than 2%. 


With this simplification the model is 


Nee WG = eit ona) (12) 
D= (N-D) (13) 
T = T(n - 1) - byT. (14) 


In state variable form where 


ee Se aa 


Xe eK oe ea) (15) 


Ke Ye Nee a yeXee (16) 
and the variable x, is defined as 


ee igi /(K, eee ele ee) 
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For a numerical example, the known constants are as- 


sumed to be A = 0.4 and y = 1.0) and the parame tenseteune 


identified are "a'' and "b'"'. Using these values the equa- 


tions are 


s 
I! 


0.4 (x, - X,) 


7S aaa. nee OR bx 


and the variable x, is defined in equation (17). 


computer equations using the scaled variables 


“fe 
I 


40x, 
x, = 40x, 
and 
Xx, 
are then 
X, = 8x, - 0.4%, 
eo = 20 = bse 


2 3 2 


where the variable x is defined by 


a 2X, 
IRE = ee, 


(18) 


(19) 


The analog 


(20) 


(21) 


(22) 


Figure 2 shows the implementation of the above equa- 


mrons and the settings of the potentionmeters for this 


scaling are shown in Table I. 


ike 
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no t 
POO] 
P002 
PO003 
PO004 
POOS 
P006 
P010 
FOLL 


POTENT IOMETERS SETTINGS 


Value 


Ore0 


Ue SW. 
0.500 
Oe 10. 


0.400 


TAB TE 
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Pot 

POULZ 
PO15 
P017 
POZL 
POS] 
PO55 


au nee 








ITI. REGRESSION ANALYSIS 


Chad 


Several attempts were made in this study to achieve 
meme desired goals for the identification system as Stated 
mioection I. Arrising in all of these attempts was the 
muoplem stated by R. C. K. Lee, [14], and illustrated in 
feenre 3, which is related to the closed-loop formed by 
the estimator and parameter identification system. This 
meonlem 1s difficult because the overall system 1S noniin- 
meee ven for linear plants. Because convergence of Chis 
Moon was not proved, and because estimation and identifica- 
tion of parameters for nonlinear cases depends strongly on 
Ge initial value*, it was decided to approach the problem 
mmcne tollowing alternative way. 

imessySteme1s noOnlimear.as detined by equatrons. (9) 
mm {11). However, the simplification performed on the 
model that leads to equations (18) and (19) yields a second- 
wader linear model with the variable x, representing an in- 
Pie. ocince it is possible to measure N and 7, which are 
mepresented in the simplified model by the variable x, and 
Mmmectace Variable xX,, mespectively, then the nom-measurable 
eeeee 1), represented by X,, Can be easily estimated. In 
@eder to make the simulation realistic, the variable x, 1s 
treated as a non-deterministic forcing function. This func- 
elon 1S generated by applying equation (22) and the resulting 

2 

PiEpis possible Lorvthiis Kind Of Ssysteiyte. reaen a 


steady-state value different from the actual value (this is 
known as a state of nonlinear stability). 
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value is corrupted with some Gaussian noise of mean zero 
obtained from the low-frequency Gaussian noise generator 
mothe Comcor-5000 analog computer. 

This simplification has the advantage that 1t opens 
the estimator-identification loop, Figure 4, because the 
mime estimated state (x,) 1S neither a tunction of the un- 
known parameter "a'', nor a function of the unknown parameter 
aD... 

Phe technique MuSed Tnmtals ws tudyven knowneasestatis teal 
/meeression analysis, consists basically of finding the best 
meee to a set of data using a least-squares criterion. It 
has been shown, [4] and [10], how the steepest descent 
method can be used on the analog computer for parameter op- 
timization and model building, and therefore this method 
Meas applied in the derivation of the regression equations 
m7) and (30). 

iteeOnGer tO personm fhe Teyressien analysis ert 1s pe: 
mmered to define S, a cost function. The cost function, 
mor simplicity, is usually chosen to be the integral of 
the square of the error, which is defined by a relationship 
that contains the parameter to be identified. The second 
meep 1s to find the rate of change of S with respect to a 
new independent variable’, t, and to minimize the cost func- 
tion using the method of steepest descent. The details of 


these methods are developed in the following sections. 


The variable t represents computer time and it changes 
more rapidly than real time. 


25 
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A. IDENTIFICATION PROCEDURE FOR THE UNKNOWN PARAMETER B 


From equation) (14) an erroresunet ion eye an be defined 


as 


e, §T - y(N-1-bT) | (23) 


ema a cost function Si als 


e 


S, ‘Gg 


1 
e? dt = f thtori-bn dt (24) 
0 Q ; 


Deis an estimate of the parameter b. 


ice pateralederavative GH Wath aweSpeet CO juamcammss 








aS, ; Jon sb 
OT ab oT 


Using the method of steepest descent applied directly 


me the objective function it is required, [4], to make 


ab JSa 
Paes | ere 
OT ab 


arse ensures knat as tT inereases tae eStimate b changes 
in such a way that S, decreases. When the estimate b makes 
3b/3t = 0 the minimum of the cost function has been reached. 

Taking the partial derivative of equation (24) with 


Aw 


mespect to b gaves 


t 
JS, e 
te 2 | WS Oso Islet} iay- it hie 
ob 


0 


and, therefore 
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a lé 1s ce 


* | 
2° = -2y J teat-y frat s rf Tdttbd rt To Gail teen) 
fale ~ 0 


0 0 0 


metting 
ee D 
eens (26) 
x, =N 


emda substituting (26) into equation (25) gives the final ex- 


pression 


AE 1 ic € 
a0 = 2 ii Kaxadt-y [ xexsdtey [ xedteb f xia : (27) 
* fo 0 0 0 


Since xX, and x, are both measurable and continuously 
patlabilwe, the identification of the unknown parameter "b"™ 
eoud also be done continuously and the input to the iden- 
tification system does not require any digital processing. 
imestO a lack of analog multipliers, however, all multipli- 
eeeons were performed in the digital computer. Figure 5 
Smews the implementation diagram for equation (27). In 
Section III-D a modification to the procedure is introduced 
muporder to avoid the need for reinitialization; this modi- 
imeation uses the properties of exponentially-mapped-past 


amerages. 


BY IDENTIFICATION PROCEDURE FOR THE UNKNOWN PARAMETER A 
ihe parameter “a can be identitied using equation (le 


Since Dis a non-measurable state and it is needed in the 
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megression equation, a Kalman filter {sas sdescriocad sino e. 
tion III-C, was used to obtain the estimate of the variable 
D.. 


Taking equa tacne( hz) 
i. D/ (K,-aT-1) 
ma defining an error criterion €, as 


- 


Eo Det N(K,-aT-1), (28) 


Teads to the cost function 


ie 
2 2 
eo [ oi dt. (29) 
0 
feos an eStimate of the parametex "a'’. The steepest descent 


algorithm can be applied to obtain the differential equation 
that defines the regression analysis for the unknown parame- 
or a". 


Using (26) gives 


E. Sua Cn roa i) ae 
and 
aS, 3 F A 
—x = — | eee (K.-a a iL) |e ite 
da da 


0 
ie 


2f [x, + X, ce = 4 na = I IgCeec grease Ghee 
0 


Therefore, 


ads 
Alm > 
tl 
Le] 
2 ¢ 
Gu 
4 
a 
to 
Pa 
Cu 
Gr 
fe 
roo 4 
(aE 
aN 
0) 
4 
tN 
ra 
w 
ar 
Gu 


De) 








iE 
2 & m. x OM = f se xX, dt |. (30) 
0 


Figure 6 shows the analog implementation of equation (30). 

As in the identificataoneotsparamcter. > a see mnem 
III-A), this implementation was augmented by using the 
properties of exponentially-mapped-past variables (see 
Section III-D). This was done in order to avoid the trouble- 


some iterative reinitialization. 


ee KALMAN FILTER 

As stated in Section II, the linear equations (15) and 
meoje are used in the Kalman filter for the estimation of the 
Beate Variable D. These equations are repeated here for 


convenience. 


= A (Ocoee (tS) 


~ 
{ 


eee Yl, = b= bX.) . (16) 


fheene application of the Kalman filter the variable x, is 
Mecated as a forcing function. 
Figure 7 shows the block diagram of the estimator for 


miessystem defined by the difference equations 


x(k=1) 


2g 


i> 


(k) + I w (Kk) 
z(k) = Hx (k) + ¥ (k) 


where k represents a point in the discrete time domain, x(k) 
meethe vector of the state variables at time ky WK) ot Saeelic 


system input which in this case is the sum of two components, 


Jie 
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the variable x, (a deterministic anpu teams Gates eae ce 
z(k) is the vector ot measured) states) a0 limes ands ae) 


is the measurement noise. The matrices ®9 


oo SD 


IT gmc Hare’, 
mespectively, the state transition matrax, thesduceribuiuen 
matrix and the measurement matrix. 

Recurrence relations for the optimal filter, [14], 


mis), are 


x(K|k) = x(k|k-1)+G(k) [g(k) -H x(k|k-1)] 


x(k+1[k) =  x(k[k), 


where x(k|k) is the estimated state vector at time k using 
observations up to and including time k, x(k]k-1) is the 
predicted state vector at time k, using observations up to 
time k-1 and G(k) tS the gain Matrix. 


The gain matrix G(k) is given by 


Gk) = P(k[K-1)H' [H PB(k|K-1)H) + RP 


P(k|k-1) is the covariance matrix of prediction error defined 


by 
Pck|k-1) © B{[x(k) - x(k]k-1)] [x(k) - x(k]k-1)]"}. 
The auxiliary equations that determine the propagation 
Sie are 
Pc(k|k) = (L- G(k) H) P(k[k-1) 
P(k+1[k) = @ P(k[k) @ + Q 


where P(k|k), called the covariance matrix of estimation error, 


is defined by 
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P(klk) 4 EC (x(k) - xCk]K)]Ex(k) - xCk|k)]7}. 


hee the covariance matrix of mthe Measmmemem Gane ic camel emis 
the covariance matrix of the system input noise, are de- 


fined as 


R = Ely(k) vy! (k)] 


Q = FElwCk) w7(k) Iz. 


The matrix form of the difference equations, that cor- 


meepond to the differential equations (15) and (16), used 


for the filter, is 


ae 2 0 008 0 





Sar x(k) + 
0 990 0 07 
0 
H = 
Ih 


fOr a Sampling interval of .02 sec. and a value of .5 for 
me parameter "b"'. 

since "b' is required for the © matrix then the iden- 
jrication of "“a'’ is done after the identification of b has 
been completed. 

The signal-to-noise ratio is specified as the ratio of 
two exponentially-mapped-past averages, the square of the 
Signal and the square of the noise signal. This average 


function is explained in the next section. 
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The computer listing of the estimator algorithm con- 
taining the complete set of equations that characgtenize the 


estimator is Within the subroutine ESTIM. 


D. THE EXPONENT TALE Y-eAPPED-PAST AVERAGE FUNCTION 

The open-loop approach wsed for the implementacioneor 
the Kalman estimator produced system equations which do not 
depend explicitly on the unknown parameter "a''. This imple- 
mentation of the Kalman estimator allows use of the regres- 
Sion analysis technique. However, this procedure allows 
meeevariable x, to reach the identifier without being fil- 
menead, Another problem still present is the iterative re- 
mitcaialization. 

To solve these problems an exponentially-mapped-past 
(EMP) average function, [13], was found very useful. The 
Eee situnction also was found to be very suitable for acting 
Moun as a filter and eliminating the overloading of the in- 
meprators used in the identification. This exponential 
weighting function decreases with early values of time and 
its implementation allows for relatively fast control of 
mEecetime Constant. It gives the advantages of easy com- 
putation and simple implementation. 

the EMPwaverage function 1S definedras 


foe) 


f(t) = af i) SH a ee ae cle 


=- 0& 


farere ult) 1S a unit step function. “This convolution imn- 


tegral of f(t) and an impulse response of a low-pass RC 
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filter, with leakage time constant RC = 1/a is shown im- 
milemented in Figitress. 

Figures 9 and 10 illustrate the changes made in order 
mooinclude the EMP tunctven anto the Simulations ihe 
Changes are the addition of input potentiometers with set- 


mimes of a. 


OS 








R=1/a 
eit) f(t) 


Figunyewoen ExXpoOncnemally-mapped-past average 
Gr Relll Ge. 
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LV pe EXECU EO Wage te 


One of the most amportante tasks prior temorrarmene 
data from the analog simulation is the verification of the 
aeerational performance of the analogs” computer. » ihas is 
Semccially critical where the analog computer used in the 
Simulation is large and has many other users, [20]. This 
task is very difficult and no attempt was made here to cover 
the whole system; however, it was found worthwhile to check 
out all components used in this simulation. 

The executive routine MAIN can be divided into two 
Emiesections: the first one performs the checkout and the 
feeona one controls the actual identification. 

The test performed by the checkout routines are: 

Paleo. - ims  SUbpnocran tests the Webyary -sub- 
ieeimnes SETPOT, SCAN and DAC. To accomplish this TESTO 
sets potentiometers and compares the desired setting with 
@iemactual setting, allowing for the inherent uncertainty 
Present in the fourth digit. In the same way the digital- 
Memanalog converters (DAC's) are tested for proper setting 
and also for any interaction among them. Any detected 
Sieeot Will cause a message to be printed and a PAUSE to be 
Gme cuted. 

Zee LEST. — Wiets: subsorcman Ces tsmche mamgeal con- 
dition setting on the plant and the analog-digital conver- 
ters and digital-analog converters. It takes 1000 samples 
of each set of converted values and prints the mean errors 


Haceone error Vvabrtances of beth tests. 
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3s TEST2. - This subroutine’ tests, themintemmrrons 
ised in the simulated nuclear. réd@eer s,s tem aeons 
ue the percentage error of Lhe Meme corate. imtce wane 
males. 

a= TESTS. -S9Nis stbyrogram tests tChemiiitral condi 
meons on the simulated plant, which allows verification of 
mime divider M003. It also runs the plant to compare the 
final value with a precomputed numerical solution. This 
M=mEGutine prints the percentage errer for beth tests. 

Se Jeo la Wits SUbroMmine periorms a test of Doth 
momettication syStems using a precomputed numerical solu- 
buon for the steady-state and prints the percentage errors 
of both systems compared with the parameter value used in 
moe numerical solution. 

The identification of the unknown parameters is per- 
memmed Dy Subroutines RUN] and RUNZ. 

ie RUND = NeGbivates the identifieation system for 
mien parameter “b'’, which is performed first because the 
Kalman filter requires the value of the unknown parameter 
Sees Ihe intérnal subroutine SAMP1, activated by interrupt 
Soeeperforms the sampling of x, and x,, Computes the deriva- 
tive of x, by a difference approximation, sets the values 
@aene reactivity input and the required products in the 
Simulated nuclear plant and identification system respec- 
mively. RUNI prints on the teletypewriter the value of “b"'. 

Zoe RUNZ e =] ACtivates the s1dewel mueacilonec)c tenmasens 


the parameter "a"'. HANDY, an internal subroutine activated 
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by interrupt 53, samples «> and) x7 Gal ls) Eo liit ees 
tilter) to obtain an estimate of 2Weepertorms sce oumbmcd 
Mmaitiplications and sets these values in the anal@pmconpucen 
for the identification system. 

The program was written primarily in FORTRAN IV for 
the XDS-9300; some assembler language statements were used 
mimorder to be able to use the library subroutine DELAY. 

A sample of the output taken under conditions of 
mieluce of the library hybrid executive routine is shown 


mea part of the Computer Output. 
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V. SIMULATION RESULTS 


Several tests were performed for each identification 
Suetem. The first test of €achesystem is eshown ine rae- 
ure 11 and Figure 12. These figures show the response of 
the identification system to a constant input whose value 
mas numerically precomputed and corresponds to a stéady- 
meate value of the plant, i.e., the digital multipliers 
Saewn in Figures 9 and 10 were furnished with these steady- 
State values. Since these runs were performed without state 
variable trajectories they do not reflect the previously 
epecitied desired objective of an identification system. 
Nometheless, for this example the identification system 
converges to the actual parameter value if supplied with a 
memeseless discrete sample representing the steady-state 
femee Of the plant variables. The time scales are un- 
labeled since the convergence time can be altered arbi- 
trarily in any noiseless situation by adjusting the system 
fain, without changing the final result. Figure 15 shows 
meemevpe Of plant transient condition used, a truncated 
map. For the identification of the unknown parameter "b" 
a low gain was used and for the parameter "a'' a high gain 
was used. These choices of low and high gain were made 
arbitrarily because for noiseless runs the procedure was 
Meersensitive to gain and e1ther high or low gains produced 
time Same satisfactory results. Here again, for the reason 
stated previously, there is no need for labeling the time 


Seares; 
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Figures 16 thru 235 show the bellavicr of (the wsyseem 
under simulated neal Jeondit nonce 

The system response in the identification of the parame- 
mer "'b'’ with two ditterent Nevers or meacmrement ance macmm 
shown in Figure 16, where the noise was added to the vari- 
feebe that, represents) the plant temperatune. x — 

Figures 7 and 8 show the system mesponse Taeeie 12¢dens 
MmmerGaliOn Of parameter “a with two daitrerent levels or 
measurement noise. Unbiased noise was added to the vari- 


Meme that represents the plant temperature, x and to the 


>? 
vetable that represents the meutron density, x,. 

Figures 19 thru 21 were obtained with biased measure- 
ment noise, and Figures 22 and 23 show typical inputs to 
fame identification system. 


Talbe ll sttmmarizes the mumerical results tor Varicus 


mms with different values of the parameters “a'' and "b". 
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TABLE [I 


NUMERICAL RESULTS 


PARAMETER IDENTIFICATION 
a b A b 
0.4 Ws 0.40 0.29 
0.4 0.4 0.40 0.39 
0.4 0.5 0.40 0.50 
0.4 mS 0.41 0.61 
0.4 0.7 Qed We 
0.4 0.8 0.39 0.81 
0.2 0.5 nai 0.50 
ies 0.5 . 0.29 0.50 
0.5 0.5 0.50 0.49 
0.6 0.5 0.59 0.50 
7 0.5 0.69 0.50 


The signal-to-noise ratios used in all these runs was 20:1. 
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Peoure 11. Identification of parameter "b"' under steady- 
state conditions. 
(2. | =e See Fas 
= te a pa ==SSS= 
ee ee ee = S| 
——— 
|= 
=== 
Saaee eee | 
= 25 
0.0- = SEES 2S=S=S=S5 
Seo 2 2a eS eas S> = == == === = 
TEMS 
Wmeepeoe 17, Identification of parameter “a'' under steady- 


state conditions. 
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mere 13. %Identification of parameter "b" under transient 
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mieure 14. Identification of parameter “a" under transient 


GOnditions. 
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Bugure 15. Reactivity function used in the simulation. 
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Figure 16. Identification of parameter "b" under noisy 


conditions. 
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meoure 1/7. Identification of parameter “a” with measurement 
noise on x, 
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Peure 18-2 Identification of parameter “a™ with measurement 
MS LSE 19) he Gey 


20 











‘ » . , } , 
2 2 { t ; ' 4 

‘ ' i : } 

BEC er ere ie 

t 

ett a — ie i 

’ ¢ ’ { .- - - t 
{ 2 4 { ZI t t { 
4 t + i i 





nel ee a a 
0 10 15 
TIME 
0.04 ae aaa ae ee 
bom pos i. 
: | | Sas Seas 
2S) Wie ae a. Ba See ae 
| | | ae ! = | | 
Petsig 7 pe | ia | | bls ? 
= eee Seat ae ste eater 
ed rn eat, aes aon eth il oo Bape geaeeronl 


; 





oa > 


pezreccenseees TT 


‘- oT 
t t i t 
| ae 
rae ns I | | Pot et EE 
Jen her } } 
(ie | on i | | 
a: | -_ i , a ee 
e 


Scere | 
= oe a 
Pee | ee eae ee 

+ Ss a {i 
PEER ee ere 









0.0m TIME Se 
Ey elo Za eee : 
Coo ae eaae i 
x ei EE Hs ia 
2 a i i | ] 
; Te i : all 
i, 0 i VOT ST or Bt 
; | ee pee 14 tbo pa cS a aie a e ge eet: oat 4 
== ; | | 
BSH HEF EE EE 
BREE E RCE | | = 
LL Shihan See oe | | | 


Figure 19. Identification of parameter ''b" under a biased 
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Expreme Diased  conditnron for parameter ‘'a'’. 
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VI. CONCLUSIONS 


From the results of thewnumerpreal cxane hemusicd sarc 
following conclusions are reached regarding the desired 
goals for the identification procedure. 

A typical response of the identification system where 
Mme initial valwe of the parameters are set to zero is shown 
i” Figures 16 thru Zl. The results indicate that the pro- 
cedure does not require a guessed value to be effective. 
Other runs were made with various initial parameter values 
and if the guessed value was close to the actual value, the 
convergence time was greatly reduced as might be expected. 
In any case, whether or not a non zero initial value was 
feea. the system did not require iterative reinitialization. 

mie Ssimplic#ey of the tdemtification procedure imple 
mentation, together with the small number of required analog 
Components, provides for a very reliable system, since ana- 
log integrators and analog multipliers are very reliable 
components. 

Mies average Convercence Sine iE wdeesa Ona -<o-louse 
ratio of 10:1 was found to be approximately 16 seconds, and 
about 10 seconds for noiseless measurements. If the signal- 
toO-noise ratio is small the convergence time is large. 

hon the exanple Gonsfeered In this study = tiesprocealnc 
identified the unknown parameters within 3% of their actual 
values with signal-to-noise ratios as low as 10:1 in the 


measured states. 


26 








The identification system for parameters b>  wasever, 
insensitive to the noise, except for biased noise (with a 
mean of 10% of the stead-state value of the state variable) ,; 
as shown in Figure 20. The extreme biased condition was 
taken as the maximum noise for which, even when the system 
mees mot actually converge, the identifier output oscillates 
within 4% of the actual parameter value. 

The identification of parameter 'a'"' also showed very 
good performance but it was more affected by measurement 
Mose On X3 than on Xo. This was expected since X3; is not 
filtered by the Kalman filter and x3 is an important factor 
mmeevetrTy term of equation (30). This problem can be elim- 


maced by using filtering in the measuring circuits for x3. 
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APPENDIX A: DYNAMIC EQUATIONS SIMPLIFICATION 


As seen in the case of an heterogeneous reactor made 
up of n independent media, the state variables are related 


By the following non linear differential equations: 


g 6 

d 

Te = @ - x % TF - e n + Dos Gals 
= 11> 

|| 

“dt oe Qe 1? AC; CliaZ) 


dT. \ 
ec; —ai- fa : 7 Dy 954075 > Ty (1.3) 


O 


The initial conditions for this set of equations are 


n(0) = a 
25) 
C; (0) fs fae - LEN. are 
T.(0) = 0. 
ae 


Por etc nUGpOse Of —tiis StTUGY tne GoOMplexity Of) tie 
model is reduced to a single time constant approximation, 


ice, 1 = 1 and j = 1, so that 


Se = (¢ - oT) ge - gent AC (2.1) 
dew 6 
ste ye ee 2 
E st = ula > 1) Sao C259) 
O 








with the following initial conditions 


A further 


morn, L.e., 


dn 
dt 
or 
1 dn 
i ale 


p., (0) 
n(0) 


— 
ry Oe 


— 
| TO 


— 
we] TO 


es 
ae 

I 
TW 


In the same way 


dc 


dt” 





OD 
WRI] O 
i 


i 
OD 
le 
1 
WI 2 


MIL 


C(0) 


T (0) 


SS 


manipulation will 


1} 
) 





preduce 2a More tecmientenn 


(ok) 


(3.2) 


(329) 








Now if 


A on A C Ap 
Nios De ae ea 
O O e B 
A 
yFtel a 8 b ee , 
then 
ae caer 
Ue (K-a yo seb Aas) 
D= AX (N - D) a2 
T = y (N-1-bT), | Cs) 
and 
N(O) = 1 TC). = 0 
D(O0) = 1 K (0) = 0. 
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APPENDIX B: REACTOR STABILITY 


As it was shown by Liapunov, “Problems general de la 
stabilite du mouvement," the general problem of stability 
of motion relative to certain quantities can be reduced in 
many cases to a solution of problems on stability of the 
mmumyial solution of the system. In theorem 1, “Direct 
Method," he stated: "If the system of differential equa- 
tions that defines the system of differential equations 
that defines the system is such that it is possible to find 
emeesitive definite function V, the derivative of which, 
With respect to time, calculated by virtue of the system 
equations satisfies the inequality dV/dt < 0, then the sys- 
tem is stable." 

iemnnonlchirs  eOmIivcS PI Oate whe tier Or eNO cic 
fuedidonary state is stable with respect to initial distur- 
bances. 

The operation of the reactor (mathematical model) is 
tmmimed by equations (3.1), (3.2) and (3.3), assuming that 
muemexternal control is disconnected and neglecting the 
Peasoening effects (since they represent a very slow phenom- 
enon) the reactivity can be considered as some function of 


meamocrature, 1.€.:; 


6K = - £(t). 


First of all it should be noted that n and C are always 
positive (and so are N and D), since they cannot be negative 


for physical reasons. Now if the initial conditions are 
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n = n(t) 
C = C(t.) 
T = TCE). 


faen the solutions 


fiagtcial conditions 


n(t) > 0 


Cc 0 


n(t) and C(t) which correspond to these 


satisfies the inequalities 


SO the stability will be investigated in the domain D{n > 0, 


mew), inherent stability in the presence of a Single group 


of delayed neutrons 


os -£(T) n - at ic (5.1) 
d 

<a = sen - AC GoD 
dt _ 5 

ee eee ae) 


which obviously admit the 


State 
T = T, some 
nen =] *+ 
O 
C= C.. 
O 


NOW eee Oe be Us 


special solution for the stationary 


LeU eS yiihe Sie ees fo. 2) 
bT (6.2) 
Cored) 


1l+b =0 are verified then 


this solution for the system is asymptotically stable; for 


physical reasons the first four inequalities of the statement 
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are always fulfilled and the fifth is a consequence of the 
fourth one. 
Considering the Liapunov function 


C 


AL n 
V - f ear + f ES) GN ar [te Gew(G7)) 
n C 
hs Ny C 


O 

miewhich the constants T,, i and C, are given by formulas 
ftorel), (6.2) and (6.3), the function V becomes zero for 

~ = Wee ie and C = Co and under the conditions given in 
Theorem (1) by M. V. Popov in [5] it is positive in the 
Gomain n> 0, C > QO. 
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Since dV/dt is negative semi-definite, it can be con- 
cluded that the solution is stable in the sense of Liapu- 
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COMPUTER OUTPUT 


The computer output shown here correspondss to duster cne 
runs where failures in the XDS-9500 library hybrid executive 


package were detected. 


Bemok IN POOG 
Beeck IN POle 
eee eR IN PO57 


More RACY YGN 
INTERACTISN 
Par ERACTI@N 
INTERACTIEN 
INTERACTIEN 
INTERACTIAN 
INTERACTIPEN 
INTERACTIBN 
INTERACTIEN 
Mey ERACTICN 
le EG RACTIEN 


BETWEEN 
RETWEEN 
BETWEEN 
BETWEEN 
3ETWEEN 
BETWEEN 
BETWEEN 
QETWEEN 
BETWEEN 
BETWEEN 
RETWEEN 


©20 
ns. 
Je) 


erz i 
T4ec 
T423 
T4eC4 
T425 
T4226 
T4e7 
T430 
(Ese 
ese 
h433 


AND 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
AND 
ANC 


T420 
T420 
peo 
T420 
T420 
T4200 
T420 
T4CO 
T420 
T4290 
T4290 


Tes } 
A/DsD/A TRUNK LINES 
PERCENTAGE ERRORS: 


MOVs SieNe: 
VARIANCES 


#202460 


©Q00000 ©00000 


een 2 
MmeGRATSRS RATE 
PERCENTAGE ERRORS: 


tated. ) 1200012 
Rave } 000232 
\Jesigs 
INT TIAL CE8NDITISNS PERCENTAGE ERRSRS: 


hEGTREN CENCENTRATISN 205493 
DELAYED NEUTRON CONCENTRATION 2000037 


65 








FINAL VALUE PERCENTAGE ERRORS: 


BEEAYED NEUTRONS CE8NCENTRATION ©44188 
TEMPERATURE °28460 
ReUTRONS CONCENTRATION Heo LoS 
ec! 4 
Pemex 7iFICATIGN PERCENTAGE ERRORS? 
PARAMETER 8 10+00000 
PARAMETER A 2950000 


66 








GN3 
(THT/19839Z2 NI Phd 18di %Xct* 
/itSW te NI TT sea AVISG fanli WSU NI SOVWNV1‘SXeT* 
4, £4459 HDLVd WAG SNNY NI SI XDEID LYHL ASZTYSA, OXSSTHTILVReO8d 2 


ZNN& WVIV3 

W YSLSKVEVd NMONHNN ASTLNSOI 

TINNY WVW1V9 

G6 YSLaghVeWd NMONANN ASTILNSC!I 

1 Sean? 2 

WILSAS NOILVIISZILNSOI LS3l 

EISat 4% 

INIA INIA ONY SNSIJDJIGNGD WILINI LIN¥4d LS3L 
elS3l V1V3 

SL1VH SHSLVSSSLNI LS3lL 

TLS31 VIV9 

NSISHYSANGD V/O"G/¥ ONV S3*D*] SYBLVYDSINI LS3L 
OL1S3i VIVO 

IV GNVY ONILLI3SS198d LS3L 

aSnvd 

(2°9)3SLIUM 

spMSLNIYd SNIV NI SNGILINYISNI GQV3N AC TOT)LNdLNAS 
JSISYNS 171V)2 

WYYDGNd NIVW 


UU OU U OU 


UO 


) OU WU U 


67 





COO0ee toy isivo 

eTft=!l O€ QC 

ANNTILNGD 

T+tN=N 

le) oy Ca ee eee eer 

SOC bee? / G2" aAxK) Say =v 

OP Gl Cole OCG ee dal* (ee kAx ea) i 
(AX fC TY IVOTIINVIS VIVO 

(O00e 11a e lay 

et‘hel 01 30 


SI DVOING@Aa a2] erica 


| SANTLNGO 
(2°GstxGlrv’s NI MGs, “XG) LVWaas 
¥ (Cl ted eo) sass 

T+NeNn 

©OOTHEZ*/S(2*eAX)SEVEV 

T OL SOCEOOO*® #37" (2*"AX)S8V) 4! 
(AX? (T)LGdIYNVIS VIVO 

(O00 "(1 jlse laa Seta 

9T*T=TI T &Q 


ONILL2S SHY3L3SWSILN3L19d 1S3L 


L3S1@q V1V9 

O=N 
LGOdHH fSSOdHth * 
TESOdMH ST SOdHt ST 2O0dHh (LTOdHY (STOdHh SeTOdHh ST TOdHH “OTOdHtY* 
S0O0dHH “GOOdHH “ HOOdHY “E0DGHY S2OOdHH*TOOCHH/IEdI VIVO 
/ECGLHe (cCEh LHe (LERLHY OGG lHe  Zewlht *9en lla 

(Gehl hchlhn (Edn lHh (SdelHe Pdr ly (OgblHy/OVd!I VLVO 
(9T)LE8dI NOISNAWIC 

(2T)IVAI NGISN3WIC 

O1S3l1 SNILNSYENS 


68 


OU UO 








ONJ 

NUN L3y 

3SNvd(O# LON) gd! 

;OLSS! NI] Seseea a (ot wei Oo 
SANITLN@S 

(OOo0od*O"lIDVG WVvO 

SANILNGD 

(ty, ONV 1 SHV. NSSML3S8 NOTLOVESLNI, OXG) LyRKuaSs 
(T)DVOL “(PF OVI (Te"9) 3L1 4m 

T#NENn 

O2 OL SO(EO00**31* (2*-AX) SEV) AI 

OAK Glo TONY oo ley oD 

Od ene 00 ( eu get) 3] 

2t‘Te O2 Gd 

(OOTO*O"T)9V0 11V5 

2t’tsI C2 8a 


SSNI7 3VG ONGWV NeILOvesSLNI sa. 


Se 


O¢ 
Ve 


UO YU 


69 








UN 
Nemes 
(S//G*ST Ad XG/S i SAINVI VAI SXSS/SS*ST® 
JS2°xG//iSHOuys SOVINSDYS di XS7 iLS3L SANIT ANNYL V/OFO/V i FXG) LVKueS 2 
SOM Sx CNG Kn re oe 
gee *QOOO0T/2X) = *O00T/2eXeSexidee( *OCOT/STX) = *O00T/2TX=SIX 
(WEX+9 OT) #9OTEWEXE CWT X= * OT) POTEnITX 
"O00 T/2KznKeX! SQOOTSTEX=ENIX 
SNNTLNGD J 
S2X#SIX+ZEXSZIX(STX*ST Xtc TXScIX 
SEX+2X tax fSTX+TX=TX 
SE2X**COTES2XfSTX**OOTHASTX 
WOOO "OCC toy Garin 
*J*! TOOV SHL iSsi TOSL ANV Teel SANTI 
NGSISYSANGD V/G9Q/¥ VIVNSAG SHL 1931 OOS S3NIND 
OPT 7 Sekesext U4 Ss bX esl x 
(SEX TS IX*O) ACY Vivo 
ABS! RARE 
OO00COTS eal 
(ODO 1 e000 Renal.) 2c alany 
COOT‘Te=i | &9O 


baS SY Gl WSISAS Sra seony so Veet shia 
SAV 14a 1SCdQesSoONvn® OCOGd Sn1q 


Be CeO on sveIN (LSS la 

O° O02¢dcXedTX=cX=TX 

(OO00 st | (O00C® (Ol Oeuee 46) 2 Garay. 
(,TiS3L s/ILVWYOS € 

Ge oo alee 

(O002*O* 20d "S° 0 “7 00da i i Garacr iaiv2 


NGISHYSANQD V/O4%GQ/¥ CNY Ss*d*!] SYyGLVYOSLNI LSsi 


PLS3L ANT LNEYENS 


70 


UUU UO 








GNa 


NonLag 


(Oot “e4 Oe rd 40°) (1) Sant eee 


Bighy > 


(S°OU a AS ala 4 1 in AS le 4 
eee ce Sree! 
Cy Cee lel 


ee 


eaod 


JOVSS3W INTYn ONV SUBGENUS BOVINEDYSd ALNdhKSD 


Mecca, Picket) eM 
SAIOVLIGA SHOGlVyOSILN! 


O78 


VW 2 
GyvalS 


Ba8 


6 Ol 69¢orlds (7 isa) ai 


(O° Le *e)Sanlatas 


VW 


SCN6D3S OF HGS SALNdWED 


3LAdWGE OQ 
(OG )42 3530 
(000°0 %200dE + {000° 0"SO0dH4) LEd1IS 
WOObeefetIoyvG 
(O* 1c" Oe =) cane 
AY 130 
OT= 
(O®"BSlIOt es (Oe Mgr Oe lorie le fT eaelergoefaorge rT) Ova 


V1WV¥9 
VV 
a2 
VW 
aay 
V3 
v1 

ay 


(/1;kGhe3 SOVINSZI&3di  XSG/ 1 aLVY SeaLVYHORBINIGSXG/icl1SSL I/D LYWeOS 
(Veo) aie! 
(2)1 NOISN3AWIO 


Ly ei Ly SOereN | 


JS 


21S3l 3NTineysns 


UU OU 


fl 








VEX SH Coxe EXeTPACY Vivo 
Q1GH VIVO 
JANTINGS 
AY 130 V1V9 
OOOT= Vd 
cood‘t=I ¢G 6a 
3iNdwad V1V9 
(CO °Weee) Se tata 


NGTLANSS WOTYSWAN NMONY 
SHL HLIM 3YVdWED ONV SSNIVA WNIS 3LNdWwsd 


AtX(T49) SLIEM 

(CA=T*) *°OT SALI X*¥O* S2O* 1) «OO L=x 
O*t/AzA 

O° c/X=X 

CAST IX SHIROV TIVO 


SYGUNS SOVINSIESd 3LNdW9D CNV SNOILIGNGD WILINI G¥Se 
AVI3SG VIVO 
O9000T= vaq 
OC isse SANT bss" diy. 
(O00h* =") *OOOTR Se (OCCT ae TCC 6 is) oy aaa 
SN@ILIGNG@D IWILINI Las 
(000d* “OT0eHt *OCOtT* *900dhn “000d * “So0d* 
Ht “CO0S* “HO0dH+ 000S* “EQ0dHh “0008* “200dMe SO002* “TOCOdHH) 1Gdli39S W1V2 
B6rarta2dfsyeco* T=EuelG 
EG*ea'TG/V/ NGWWSD 
SNAVA AWAITS GNV SNOTLIONOOD VWILINT W378Gudd 1S3al 


ELSa3i ANITLNGYsns 


UU U O 


UU U 


UU U 


UL: 








CN3 

(/S*ST sxe’ NOTLIVYLN3GD* 

NOD SNOMLNSNIOXS/S*STS“X6T ls SHNLVeSdW3L i TXS/S*ST AY XT NGI LVUINGIN® 
G2 SNGYLNSAN GSAV1901'XS/  tSeOsyH3 SSVLNBIYNSd SANIVA TWNIS i/)LVWHOS E 
(/S*Sts’s * 

NOILYSLNSINGD NOYLNSN GC3AV13G. !XSG/G*STs!XE $4 NOILVYLNSINGD NGULN* 
INi?XS/i12SY¥OHNys SOVINSIUSd SNGILIGNSD WILINT w/s€1S3L w/)LvWues 7 


Nenlay” 
EX 2X LAWS ers aM 


O* O00 Tx 
(EG/(EX*EG) SAVEEX! *OOT*(2ZE/ (2X228) )SEVAPZXfOOTH(TSG/(TX=T&))SEV=TxX 


EX*#*OT EEX fEXw*OTeacx i TX#*OT==TX 
O*CSEXBEX (0% 4 /EX2CXK fO% H/T XETX 


IS 








JOVSS3W LNIYd ONV YGNYS SOVINSADYSd JLNdWed 


g§ 40 SNA vay 


Chay eter” dy 2) 

SANTINGS 

KV AEG WY 

OCOT=s Vici 

OO002*T=1 Tt 8d 

3LAdWSD 11V39 

(Led ©8997 12 Ge 11 oe Ne) Oy Cena oD 
EGeZ2deLSnGDvV*TO*e=Z42 IL 
CGxCox*OT#LSNGDV*TO*?=92 11 
O*O02=G¢Z2 11 
RErZeTHeLSNGIveTO*==2€2 IL 


MSISTLN3GI SHL SLNI STVNOIS LAdNI LSS GNV 31NdhkOd 


(l* 

SNODV ET ZOdHh “LSNOQDV “LTOdHH “LSNODV “STOdHH (LSNEQVOTTOdHH) LOd1L3S W1V9 
(LSNGIV* 

CLGOdKh “LSNOOYV §SSOdH+ £0606 "0 “ESGOdHH fLSNGIV STGOdHt) LOd139S VIVO 

(OM *GulceiNina bas salaiy 2 


NOTLVII4SILNSG] LNSYSS5I10 SHL Yw8@4 YOOX SSHILIMS GS ANID 13S 


9°02 1SNQDYV 

SANTLNGO 

a1ViS AGVals apeesos Not ies 

TWOIMSWON NMONH VY ONISN @ &S3LEWVYEVG NMENANA GHL 46 NOTLVITSTLNGQI 
LSNGIV/E/ NOWWED 

ESLSE*TS/V/ NOWWGE) 

W3SLSAS NOTLVIIAILNACGI LS4l 

4#1S3i SNILneuvens 


GUO UOU U 


OO YU 


OO YU 


74 








U.S. DEPARTMENT OF COMMERCE 
National Technical Information Service 
9285 Port Royal Road 
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November 29, 1972 


United States Naval Postgraduate School 72-5 
Monterey, California 93940 


Re AD 738 874 Solution of a Nuclear Reactor Parameter Identification 
Problem ; 


Gentlemen: 


We are unable to provide the public with satisfactory copy of referenced 
document for the reason(s) checked below. Please help us by sending 
replacement copy together with this letter by AIRMAIL to address shown 
below. If you are unable to furnish replacement, please return this 
form by AIRMAIL indicating where the necessary replacement copy is 
obtainable. Your cooperation will be appreciated. 


1. These pages are missing, please furnish: Page 75. 
[ ] 2. These pages are not reproducible, please send good pages: 
a 3. Document is not reproducible, please send good copy. 


(| 4. These pages have been mutilated, please send replacements: 


[] 5. Other 


Input Section 

National Technical Information Service 

Department of Commerce 

5285 Port Royal Road 

Springfield, Virginia 22151 Tos: Mr. Thott ‘ 


I regret that our Library copy 
is also incomplete in that it lacks 


(A page 75, too. We have contacted the 
a advisor in a search for the original 
4, | THOTT but it appears that it has either 
Chief, Input Section been discarded or was retained by the 
703-321-8517 author, who is an officer in the 


Peruvian Navy and Zz has now returned 


to his country. an Koh 
Leg R. Luckett 


Librarian 
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